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Résumé 



Ces notes sont issues du cours que j'ai dispensé aux étudiants de première année ins- 
crits à la faculté des sciences de l'université Mohammed V pour le master "physique 
mathématique". Elles sont basées pour lessentiel sur les sources énumérées dans la 
bibliographie. 
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Chapitre 1 



Introduction 



Classiquement, il existe deux types de systèmes dynamiques. Il y a d'abord le mou- 
vement d'une particule ou d'un corps rigide avec un nombre fini de degrés de liberté, 
qui peut être décrit par un nombre fini de coordonnées. Et puis, il y a des systèmes 
physiques pour lesquels l'ensemble des degrés de liberté est non dénombrable. Ces 
systèmes sont décrits par des champs. Des exemples courants de champs classiques 
sont les champs électromagnétiques décrits par E (x, t ) et B (x, t) ou de façon équivalente 
par les potentiels ((p(x,t),A(x,t)). De même, le mouvement d'une corde unidimen- 
sionnelle est également décrite par un champ (p(x, t). Ainsi, alors que les coordonnées 
d'une particule ne dépendent que du temps, les champs varient continuellement en 
fonction de certaines variables de l'espace. Par conséquent, une théorie décrite par des 
champs est généralement vue comme une théorie des champs à D + 1 dimensions, où 
D représente le nombre de dimensions spatiales. Par exemple, une théorie décrivant 
les déplacements de la corde unidimensionnelle constituerait une théorie des champs 
de dimensions 1 + 1 alors que la théorie des équations de Maxwell, plus familiers, peut 
être considérée comme une théorie de champs à 3 + 1 dimensions. Dans cette perspec- 
tive, il est alors clair que la théorie décrivant le mouvement d'une particule peut être 
particulièrement considérée comme une théorie de champs à + 1 dimensions. 

L'intégrale de chemin est une généralisation à un nombre infini de variables, représenté 
par des chemins, des intégrales ordinaires. Elle vérifie les mêmes propriétés algébriques 
de celle-ci, mais présente en plus des propriétés nouvelles. 
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La puissance des méthodes de l'intégrale de chemin réside dans le fait qu'elle met, très 
explicitement, en correspondance les mécaniques classique et quantique. Les quantités 
physiques s'obtiennent en moyennant sur tous les chemins possibles. La caractéristique 
intéressante de l'intégrale de chemin est que dans la limite semi-classique h — >• les 
chemins dominant l'intégrale se trouvent dans un voisinage du chemin "classique", ce 
qui permet de bien comprendre le principe de correspondance. 

La formulation de la mécanique quantique basée sur l'intégrale de chemin est plus 
compliquée du point de vue mathématique. Cependant, elle est bien adaptée à l'étude 
de systèmes à grand nombre de degrés de liberté où un formalisme basé sur l'équation 
de Schrôdinger perd toute son utilité. Elle permet ainsi une analyse aisée de la théorie 
quantique des champs et de la mécanique statistique. 

1.1 Notes historiques 

En 1933, Dirac a observé que l'action joue un rôle central dans la mécanique classique 
(il considère que la formulation lagrangienne de la mécanique classique est plus fon- 
damentale que la formulation hamiltonienne), mais qu'il ne semblait avoir aucun rôle 
important dans la mécanique quantique. Il a ensuite avancé que cette situation pourrait 
être corrigée si le propagateur de la mécanique quantique "correspondait à" expiS/h 
où S est l'action classique évaluées le long du chemin classique. 

En 1948, Feynman a développé la suggestion de Dirac, et a réussi à dériver une troisième 
formulation de la mécanique quantique, basée sur le fait que le propagateur peut être 
écrit comme une somme sur tous les chemins possibles (pas seulement les classiques) 
reliant le point initial et le point final du mouvement d'une particule. Ainsi, alors que 
Dirac n'a considéré que le chemin classique, Feynman a montré que tous les chemins 
contribuent. Pour l'anecdote, l'article original de Feynman qui a jeté les bases du sujet 
a été rejeté par Physical Review ! 
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1.2 Notations 



Dans ce cours, nous allons considérer les théories non-relativistes et relativistes. Pour 
le cas relativiste nous allons utiliser la convention de Bjorken-Drell : les coordonnées 
contravariantes sont supposées être 

x?' = (f,x), p = 0,1,2,3, (1.1) 

tandis que les coordonnées covariantes sont 

x F = y] } i V x v = (f,-x). (1.2) 

Tout au long de ce cours, nous posons c = 1. La métrique covariante est donc diago- 
nale : 

r\ w = (+,-,-,-) (1-3) 
la métrique inverse ou contravariante a la même forme 

V ¥V = (+,-,-,-)• (1.4) 
l'élément de longueur invariant est 

x 2 = x v xV = 7] }IV X^X V = t]^x F x v = t 2 - x 2 . (1.5) 

Les gradients sont également obtenus à partir de (1.1) et (1.2) : 

3 fd 



3 « = ^=U' V| ' (1 - 6> 



dXfi \ dt' 



le D'Alembertien prend donc la forme 



□ = d% = tjWdjA, = ^ - V 2 . (1.8) 



1.3 Fonctionnelles 

Dans l'étude de n'importe quel système dynamique, on a souvent affaire avec des 
fonctions dont la "variable" est une fonction. Ces fonctions composées s'appellent des 
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fonctionnelles. Si nous considérons le mouvement d'une particule plongée dans un 
potentiel à une dimension, alors le Lagrangien est donné par 

i 

L(x,x) = -mx 2 - V(x), (1.9) 

où x(t) et x(t) sont la coordonnée et la vitesse de la particule. La fonctionnelle la plus 
simple à laquelle on peut penser en présence d'un lagrangien est l'action définie par 

S[x] = / dtL(x,x) (1.10) 
Jt, 

Notons que, contrairement à une fonction dont la valeur dépend d'un point particulier 
dans l'espace de coordonnées, la valeur de l'action dépend de la trajectoire tout le long 
de laquelle l'intégration est effectuée. Pour les différents chemins reliant le point initial 
x(ti) au point final x{tf), les valeurs de la fonctionnelle action seront différentes. 

Ainsi, une fonctionnelle a la forme générique suivante 

F[f] = [ dxF(f(x)) (1.11) 



La notion de dérivée peut être étendue au cas de fonctionnelles de façon naturelle à tra- 
vers la notion de fonctions généralisées (ou distributions). Ainsi, on définit la dérivée 
fonctionnelle ou la dérivée de Gâteaux par la fonctionnelle linéaire 



la quantité s'appelle dérivée fonctionnelle de F[f]. Elle est définie par 

^ [/(*)] = Hm F(f(x)+eS(x-y))-F(f(x)) 
ôf{y) e^o e 

Il s'ensuit de cette définition que 



ô(x-y). (1.14) 



où ô(x) est la fonction (plus exactement distribution) de Dirac. 

La dérivée fonctionnelle satisfait toutes les propriétés de la dérivée, notamment 



5f{xY ' LJi Sf(x) Sf(x) 

5 m\m\) = ô Mm]+w)Wê 
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En outre, une fonctionnelle peut être développée en série de Taylor comme suit : 

F [f] = J dxP (x) + J dx 1 dx 2 P 1 (x lr x 2 )f(x 2 ) + J dx 1 dx 2 dx 3 P 2 (x lr x 2r x 3 )f(x 2 )f(x 3 ) + . 

(1.16) 



ou 



Pl{Xi,X 2 



P 2 {xi,x 2 , x 3 ) 



Pq{x) = F(/(x))| /(t)=0 , 



f(x)=0 



*f(*2) 

1 S 2 F(f( Xl )) 



(1.17) 



2\ôf(x 2 )ôf(x 3 ) 



f(x)=Q- 



etc. 



1.3.1 Examples 

(i) Soit F(f) = f{x) n , où n est un entier positif. 

F[f]= JdyF(f(y)) = j dy/(y)« . (1.18) 



Alors, 



mf(y)} = lim F(f(y)+eô(y-x))-F(f(y)) 

Sf(x) e^O £ 

= lim (/(y)+e%-x))"-/(yr 

e^O £ 

= lim fivT + nefjyr^ôjy - x) + 0(e 2 ) - /(y) M 

e^O £ 

= nf{y) n - l ô{y-x). (1.19) 



Par conséquent, nous obtenons 



^[f\ _ f d / f [/(!/)] 



ôf(x) J J Sf(x) 

= Jdynf{ V y-H{y-x) 

= nf(x)"- 1 (1.20) 



(ii) Considérons ensuite l'action (1.10) 



S[x(t)]= [ f dt' L(x{t'),x(t')), (1.21) 
Ju 



avec 



L(x{t),x(t)) = Irtnx(t) 2 - V(x(t)) 
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T(x(t)) - V(x(t)). (1.22) 



Nous obtenons sans difficulté 



ôx{t) e 

= V'{x{t'))ô{t' - t), (1.23) 



où V'(x(t')) = dV(x(t'))/dx(t'). De même, 

gw)) = lim nx(n+ejrô(t'-t))-T(x(t')) 

ôx{t) e^o e 

= mx(t')-^ô(t'-t). (1.24) 

Il est clair maintenant que 

SL(x{t') f x{f)) ÔT{x(t') -V(x(t'))) 
ôx(t) ~ ôx(t) 

= mx(t')-^ô(t' - t) -V(x(t'))è(t' - t). (1.25) 



Nous obtenons donc pour f,- < t < tt 

ÔS[x(t)] _ dt ,ÔL{x{t'),x{t')) 



ôx{t) Jti ôx(t) 

= j\t' (mx{t')-^S(t' -t)-V\x{t'))ô(t' -£)) 

= -mx{t) - V'{x{t)) 

à 3L(*(f), x (f)) + BL{x{t),x(t)) (L26) 



dt dx(t) dx(t) 

Ce dernier membre rappelle l'équation d'Euler-Lagrange. En fait, on peut noter que 

àS[x(t)} d dL(x(t),x(t)) dL(x(t),x(t)) 
Sx(t) dt dx(t) dx(t) K ' 

n'est autre que l'équation d'Euler-Lagrange écrite comme extremum de la fonction- 
nelle action. Ceci constitue exactement le principe de la moindre action exprimée dans 
le langage des fonctionnelles. 
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1.4 Un bref rappel de mécanique quantique 



Dans cette section, nous allons décrire très brièvement les éléments essentiels de la 
mécanique quantique. L'approche conventionnelle à la mécanique quantique commence 
par la formulation hamiltonienne de la mécanique classique en considérant que les ob- 
servables sont des opérateurs qui ne commutent pas. La dynamique, dans ce cas, est 
donné par l'équation de Schrôdinger dépendante du temps 

m a M!)l = H |*(f)>, (i.2 8 ) 

où H est l'opérateur hamiltonien du système. De façon équivalente, la fonction d'onde 
d'une particule, dans le cas uni-dimensionnel, satisfait 

fo tyfcO = H (x)tp(x,t) 

h 2 d 2 \ 

+ V(x))f(x,t), (1.29) 



2m dx 2 
où l'on a posé 

tp(x,t) = (x\ip(t)), (1.30) 

\x) désignent les états de base dans la représentation de coordonnées. Cela définit alors 
l'évolution dans le temps du système. 

L'object principal de la résolution de l'équation de Schrôdinger est la détermination de 
l'opérateur d'évolution qui génère la translation de ce système dans le temps. A sa- 
voir, l'opérateur d'évolution transformant l'état quantique à un temps passé t 2 à l'état 
quantique à un temps avenir t\ : 

\ip(.h)) = U(h,t2)\ip(t 2 )). (1.31) 

Clairement, pour un hamiltonien indépendant du temps, on voit, à partir de l'équation 
de Schrôdinger 

U(ti, h) = 0(fi - t 2 )e~T' {h ~ tl)H (1.32) 
où 9 est la fonction échelon de Heaviside. 

Il est évident que l'opérateur d'évolution temporelle n'est autre que la fonction de 
Green pour l'équation de Schrôdinger vérifiant l'équation différentielle : 

~ U(h, h) = ihS(h - h). (1.33) 



(Note, la dérivée au sens des distributions de 9 est égale à la distribution de Dirac S). 
Pour déterminer l'opérateur U, il faut calculer ses éléments de matrices dans une base 
donnée. Donc, par exemple, dans la base des vecteurs propres de l'opérateur position 
X : {\x)} définie par 

X\x) = x\x), (1.34) 

nous pouvons écrire 

(xi| U(ti,t 2 ) \ X2) = U(t 1 ,x 1 ;t 2 ,x 2 ). (1.35) 

Si on sait calculer complètement la fonction U(ti,xi;t 2 ,x 2 ), alors l'évolution dans le 
temps de la fonction d'onde tp peut être exprimée par 

xp(xi,t x ) = (x 1 \U(t 1 ,t 2 )\y(t 2 )) 

= j àx 2 {x l \U{t l ,t 2 )\x 2 ){x 2 \ip{t 2 )) 

= J dx 2 U(ti,x 1 ;t 2 ,x 2 )tf(x 2 ,t 2 ). (1.36) 
Dans ce calcul, nous avons utilisé la relation de fermeture : 

J dx\x){x\ = 1. (1.37) 

Le schéma adopté dans la discussion précédente est celui de Schrôdinger, où les états 
quantiques \tp(t)) sont dépendants du temps tandis que les opérateurs ne le sont pas. 



Schéma 


Schrôdinger 


Heisenberg 


États 


dép de t 


indép de t 


Opérateurs 


indép de t 


dép de t 



D'autre part, dans le schéma de Heisenberg, où les états sont indépendants du temps, 
nous pouvons écrire par (1.31) 

\f) H = \f(t = 0)) s = \ip(t = 0)) = e* tH \ip(t)) = e* tH \ip(t)) s . (1.38) 

Dans ce schéma, l'opérateur position est lié à celui dans le schéma de Schrôdinger par 

X H (f) = er. tH Xe-T' tH . (1.39) 
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Les états propres de ces opérateurs satisfont 

X H (t)\x,t) H = x\x,t) H , (1.40) 



et par conséquent 



Pour t\ > t 2 , il est clair que 



\x,t) H = er- tH \x). (1.41) 



H{xi,h\x2,t 2 ) H = (xi\e T^ H e t^ H \x 2 ) 
= ( Xl \e-^-^ H \x 2 ) 

= (Xi\U(ti,t Z )\X2) 

= U{h,xi;t 2 ,x 2 ) (1.42) 

Cela montre que les éléments de matrice de l'opérateur d'évolution temporelle n'est 
autre que l'amplitude de transition entre les états de base dans le schéma de Heisen- 
berg. 
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Chapitre 2 

Intégrale de chemin en mécanique 
quantique 

2.1 Formules de base de la mécanique quantique 

Avant de déterminer la représentation par l'intégrale de chemin de l'amplitude de 
transition U(tf ,Xf,ti,Xi), récapitulons quelques formules de base. Considérons, pour 
simplifier, un système uni-dimensionnel. Les états propres de l'opérateur position X 
forment une base orthonormée : 

(x\x') = ô{x — x'), 

Jdx\x)(x\=l. (2.1) 

Les états propres de l'opérateur impulsion P forment également une base orthonor- 
male : 

P\p) = p\p), (2.2) 

et 

W) = Kv~v')> 

j dp |p)(p| =1. (2.3) 
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Le produit scalaire (x\p) donne les éléments de matrice de l'opérateur de passage entre 
les deux bases et (|p)}. En fait, on peut calculer 

{p\x) = -jL=e-fa = {x\p)* . (2.4) 

Ces relations définissent des transformations de Fourier. Plus précisément, à l'aide de 
la relation de fermeture, on peut écrire la transformation de Fourier d'une fonction 
comme 



f(x) = ( X \f) = y dp (x\ P ) ( P \f) 

dpe*P x f(p) 



1 



\/2nh 
1 

\/2n 



dke lkx f(k) (2.5) 



ou 



f(k) = Vhf(p) 



V2nh 
1 



'2n 

Ici, k = p/h peut être vu comme le nombre d'onde. 

Dans le schéma de Heisenberg, les états quantiques vérifient 

H (x,t\x',t) H = (x\e-i tH e* tH \x') 



J dxe-y x f(x) 
dxe~ lkx f(x). (2.6) 



(x\x') = ô(x- x'), (2.7) 



et 



J dx \x,t) H (x,t\ H = J dxeh tH \x) (x\ e h tH 

= e* tH J dx \x) (x\ e~T' tH 
= ei m \e-ï m = \. (2.8) 

Il faut noter que l' orthonormalité ainsi que la relation de fermeture ne sont valables 
dans le schéma de Heisenberg que pour des temps t égaux. 
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2.2 Ordre des opérateurs en mécanique quantique 



Les opérateurs en mécanique quantique ne commutent pas. Par conséquent, le hamil- 
tonien classique est sensé "devenir" un opérateur quantique : 

H(x,p) — >H(X,P). 

Cependant, ce passage ne spécifie pas l'ordre d'écriture de l'analogue quantique d'un 
terme comme xp par exemple. Est-ce que c'est XP ou, plutôt PX. Comme cet ordre 
n'importe pas dans la mécanique classique, il est crucial en mécanique quantique, et, 
à priori, il n'est pas clair qu'un ordre donné doit correspondre à la théorie quantique. 
c'est le problème de l'ordre des opérateurs et, malheureusement, il n'existe aucun prin- 
cipe bien défini dans ce sens. Il y a cependant quelques conventions que l'on adopte 
en général. Dans l'ordre "normal", on ordonne les produits des opérateurs X et P tels 
que les impulsions se mettent à gauche des coordonnées. Ainsi 

xp — > PX 

O.N 

px — > PX 

9 O.N j.,.9 

x z p — > PX 1 
xpx — > PX 1 

etc. Cependant, il existe une autre convention plus satisfaisante et qui est plus lar- 
gement adoptée. Il s'agit de l'ordre de Weyl qui consiste à symmétriser les produits 
d'opérateurs dans toutes les combinaisons possibles avec le même coefficient de pondération. 
Ainsi, 

xp ^ ±(XP + PX) 

px ^ i(XP + PX) 

x 2 p ^ i(X 2 P + XPX + PX 2 ) 

xpx ^ ±(X 2 P + XPX + PX 2 ) 

etc. Pour l'ordre normal, c'est un exercice facile de montrer que pour tout hamiltonien 
quantique obtenu à partir d'un hamiltonien classique H(x,p) on a 

(x'\ : H : . N |x) = J dp (x'\p) (p\ : H : . N \x) 

= j ^e-rM— r ) H (x,p). (2.9) 
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On peut de même montrer que sous l'ordre de Weyl, cette relation s'écrit 

(x'| : H : w \x) = J H (^p) ■ (2-10) 

Comme on le voit, les éléments de matrice du hamiltonien sous l'ordre de Weyl conduit 
à ce qu'on appelle la prescription du point milieu (mid-point) et c'est ce que nous allons 
utiliser dans toutes nos discussions. 



2.3 Calcul de l'amplitude par l'intégrale de chemin 

Nous sommes maintenant prêts pour le calcul de l'amplitude de transition (ou propa- 
gateur). Rappelons que dans le schéma de Heisenberg, pour tt > t ir nous avons 

U.{tf,Xf,ti,Xi) = u \Xf,tf\Xi,ti) H . 

Divisons l'intervalle entre les temps initial et final en N segments infinitésimaux de 
longueur e. Soit 

te — t{ 

e = -V- < Z11 > 
En d'autres termes, pour la simplicité, nous allons discrétiser l'intervalle de temps et, 
à la fin de nos calculs, nous passerons à la limite e — > et N —> oo. Un instant in- 
termédiaire t n peut donc s'écrire 

t n = ti + ne, n = 1,2,. . .,N - 1. (2.12) 

Utilisant la relation de fermeture (Eq. (2.8)) pour tous les instants intermédiaires, nous 
obtenons 

\l{tpxf,t i ,xï)= Um J dxi • • ■dx N - 1H (xf,tf\x N - 1 ,t N _i) H 

N^oo 

x fi (xN-l,tN-l\xN-2,tN-2)H ' ' 'H ( x l, h\*h k) H • 

Notons que tout produit scalaire dans cette formule admet la forme 

H {Xn,tn\Xn-l,tn-l) H = M <r*(*»-'*-l)H U„_l) 



(2.13) 



ix n I 6 f < \x n — \. 



2nh 
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e f,Hn^n *n-i/ 2 'P" J (214) 



Ici nous avons utilisé la prescription du point milieu de l'ordre de Weyl (2.10). 

Substituons cette forme du produit scalaire dans l'amplitude de transition, nous obte- 
nons 



U(tf,Xf}ti,Xi) = lim / dxi • • • dx^-i 

1 1 e— >0 J 



N->oo 



dpi dp N 
2nh ' ' ' 2nh 



x e 



f, Ln=l ( Vn (*n -X„-i) -eff( - + \ n 1 ,p„ 



(2.15) 



où l'on a posé xq = x\ et x^ = Xf. C'est le forme la plus grossière de l'intégrale de 
chemin de Feynman et est définie dans l'espace de phase du système. Le facteur de 
phase dans le membre de droite de la dernière équation s'écrit dans la limite e — > oo 



N 



n^oo n=l v ' v 



hJ ti 



dt (px — H(x, p)) 
dt L . (2.16) 



En d'autres termes, il est proportionnel à l'action. Pour obtenir la forme la plus fa- 





t. Ê 2E 3E 



FIGURE 2.1 - L'amplitude comme somme sur tous les chemins N-segmentés. 

milière de l'intégrale de chemin impliquant le Lagrangien dans l'espace de configura- 
tion ( c'est à dire une intégrale sur les variables x uniquement), nous considérons la 
classe des Hamiltoniens qui sont quadratiques par rapport aux impulsions p. Choisis- 
sons 

H(x,p) = ?- + V(x). (2.17) 
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Dans ce cas, nous avons par l'équation (2.15) 

dpi dp N 



U(tf, Xf-, t{, Xi) = lim J dxi ■ ■ ■ dx^-i 

N->oo 



2nh 2nh 

xe h n ^V^ e > "" v * (2.18) 



ieyN („ ( x n-*n-\ \ vl v ( x "+ x n-l 



Les intégrales sur les impulsions sont des intégrales gaussiennes que l'on peut facile- 
ment calculer. Utilisant la formule 



on trouve 



2nh 6 J 2nh 6 

2 



Remplaçant en plus dans l'équation (2.18), nous obtenons 

/ m \N/2 



U{t f ,x f ,ti,Xi) = Hm(— ) 

N— >co 

dx\ ■ ■ ■ dxjv-i e 
= Aj ' ®[x{t)]e*ti dt ^ mi2 - V{x) ) 

= A j ®[x{t)]e L n s ^\ (2.21) 

où A est une constante indépendante de la dynamique du système et S [x(t)] est l'action 
classique donnée par (1.10). C'est l'intégrale de chemin de Feynman pour l'amplitude 
de transition en mécanique quantique. 

L'intégrale de Feynman stipule que l'amplitude de transition entre un état initial et un 
état final est simplement la somme sur tous les chemins, reliant les deux points, du 
facteur exp(i/h S[x(t)]). Cette somme est traduite par l'intégrale sur les coordonnées 
spatiales avec une mesure d'intégration donnée par 3>[x(t)]. 
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Exercices. 

1. Equivalence de la méthode d'intégrale de chemin avec l'équation de Schrôdinger. 

2. Montrer que l'amplitude de transition pour une particule libre s'écrit 



U(t f ,x f ;t f ,t t )= I 2nm{ —- 



1 / \2 

e~ h 2(t f^ . (2.22) 



(a) Montrer que ce propagateur est une solution de l'équation de Schrôdinger. 
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Chapitre 3 



Intégrale de chemin pour l'oscillateur 
harmonique 



Comme exemple important d'application de l'intégrale de chemin, nous allons considérer 
l'oscillateur harmonique uni-dimensionnel qui est un système exactement soluble. Un 
oscillateur en interaction avec une source externe est décrit par le Lagrangien 

L = -mi 2 — -mco 2 x 2 + Jx. (3.1) 
2 2 

La source dépendante du temps J(t) peut être vue, par exemple, comme un champ 
électrique si l'oscillateur est supposé acquérir une charge électrique. La limite de l'os- 
cillateur harmonique libre s'obtient de ce système par passage à la limite }(t) — > 0. 
L'équation d'Euler-Lagrange appliquée à l'action (l.lO)donne la trajectoire classique : 

ce qui donne 

mx c \ + mco 2 x c \ — J = (3.3) 

La forme générale de l'amplitude de transition est, comme nous l'avons vu dans le 
chapitre 2, est donnée par 

U{t f , x f ; t f , U) = aJ ®[x{t)\ ei s ^ (3.4) 
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FIGURE 3.1 - Fluctuations autour de la trajectoire classique. 

Pour évaluer cette intégrale fonctionnelle, notons que l'action est quadratique dans les 
variables dynamiques x(t). Par conséquent, si on définit 

x(t) = x c \ + rj(t), (3.5) 

alors nous pouvons développer l'action par rapport au chemin classique comme 

SS\x] 



S[x{t)} =S[x d + t](t)} =S[x d }+ j dtrj(t 



Sx(t) 



x=x cl 



(3.6) 



Nous notons, d'après (3.2), que la trajectoire classique est un extremum pour l'action. 
Donc, nous avons 

= 0. (3.7) 



Sx(t) 

Par conséquent, nous pouvons écrire l'équation (3.6) comme 



(3.8) 



Ôx(t\)ôx{t2) 

La dérivée fonctionnelle dans cette équation peut être évaluée pour l'action de l'oscil- 
lateur harmonique. Un calcul direct faisant appel aux propriétés des dérivées fonction- 
nelles (voir chapitre 1) donne 

S[x] — S [x d ] + ï jj dt m (V - a; 2 // 2 ) (3.9) 

Les variables tj (t) représentent les fluctuations quantiques autour du chemin classique, 
c'est à dire elles mesurent la déviation d'une trajectoire par rapport à la trajectoire 
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classique. Puisque les extrémités des trajectoires sont fixés, les fluctuations satisfont 
aux conditions de bords : 

niM) = y{t f ) = o (3.io) 

Il est clair que la sommation sur tous les chemins est équivalente à celle sur toutes les 
fluctuations possibles sujettes à la condition (3.10). Par conséquent, 



U(t f ,x f ;t u x t ) = a] toie^+^k*^ 



AeisM f tojebSÏlàtmif-uAf) (3 . n) 



Nous remarquons que l'exposant dans cette intégrale est quadratique par rapport aux 
variables d'intégration. De telles intégrales peuvent être évaluées de plusieurs façons. 
Nous proposons ici une méthode basée sur un développement de Fourier. 

Dans l'intégrale (3.11), on redéfinit la variable temporelle comme 
Donc 

U(t ff Xf f t if Xi) = Ae* s ^ j ^ e iIo rdt K-»V) (3.12) 
où T = tf — t{. La variable t](t) satisfait à la condition aux limite 

ij(0) = ti(T)=0. (3.13) 

La fluctuation tj peut, par conséquent, être développée dans une base de fonctions 
trigonométriques. Par exemple 



^(0 = E fl * sin (~y~) > n = l,2,...,N-l 



(3.14) 



N étant toujours le nombre de subdivisions de la trajectoire segmentée du chapitre 
précédent. Remplaçant, on trouve 



9 / \ t— ■ f T , nmn ( nnt\ 
dtrj (t) = J Q df^-^cos^— J 



cos 



mnt 



îiïtJo T 2 \ T J V T 

^E(f) 2 ^ (3-15) 
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où nous avons utilisé la propriété d'orthogonalité de la fonction cosinus. Nous trou- 
vons de même, 

n 

Notons que l'intégrale sur les variables rj(t) est équivalente à l'intégration sur les va- 
riables a n qui sont les coefficients de Fourier . On peut par conséquent écrire 

U(tf,xf, ti, xi) 

= lim A'eft s ^i] dai ■ ■ • da N _! e 2h Ul=l ^ T 

e->0 J 
N^oo 

;„ r , r imT v-N-1 (( nn\ 2 ,,2\ n 2 

= lim A'eî^iJ da x ■ ■ ■ da N _i e ^ E »=i _a; ) a « (3.17) 

N^oo 

Ici, nous notons que le Jacobien du changement de variables de tj aux coefficients a n a 
été "absorbé" dans la constante A' dont la forme explicite sera déterminée par la suite. 

L'amplitude de transition est, dans ce cas, un produit d'intégrales de Fresnel que l'on 
peut aisément évaluer. En fait, nous avons 

-1/2 



i rN-l(Tfn7tV ma 2_T mcv 2 a 2 



œ 2 



f4nih\ 1/2 T f /o;T\ 2X ~ 



-1/2 



\—t) — 1_ — • (3 " 18) 
Substituant cette forme dans l'expression (3.17), nous obtenons 

U(t f/ x f t l/Xl ) = Um A"^M f[ ( x _ MX\ f (3 . 19) 

N->oo «=1 \ V ' J 

la constante A" "absorbe" toutes les constantes. Si nous utilisons maintenant l'identité, 

i— r ( u 2 \ sin tcu ,„ 

nous obtenons après passage à la limite N — > oo 

U(t f ,x f ;t t , Xl ) = lim A"e^ • (3-21) 
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Nous pouvons calculer la constante A" simplement par noter que lorsque co — >• , l'os- 
cillateur harmonique se réduit à une particule libre. En fait, identifiant l'équation (2.22) 
et l'équation (3.21) pour co = 0, nous obtenons 



lim A" = f^-V. (3.22) 
\2mhTJ v ' 

N->oo 

La forme complète de l'amplitude de transition peut finalement s'écrire sous la forme 



i 

>2nih sino;7\ 

Quant à l'action classique S[x c i] de l'oscillateur harmonique, elle peut être directement 
déduite du Lagrangien : 

S fol] = rp { (tf + 4) ^ ~ 2X l X f} t 3 - 24 ) 
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Chapitre 4 



Méthode de l'intégrale de chemin pour 
la théorie quantique des champs 

4.1 Introduction 

Nous allons voir dans ce chapitre l'application la plus intéressante des intégrales de 
chemin. Comme nous l'avons mentionné dans l'introduction, cette méthode est surtout 
bien adaptée à l'étude de systèmes à très grand nombre de degrés de liberté, voire une 
infinité comme pour la théorie des champs. 

Notre objectif est de présenter une explication élémentaire des mathématiques derrière 
les intégrales de chemin de Feynman appliquées à la théorie quantique des champs 
(TQC). Notre prototype est une TQC connue sous le nom de théorie des champs sca- 
laires avec terme d'interaction en Quoique l'exemple est élémentaire, il fournit suf- 
fisamment d'informations pour comprendre l'application des intégrales de chemin aux 
systèmes de plusieurs particules en interaction. 

Le passage de la mécanique quantique à la TQC est direct, mais le concept sous-jacent 
est un peu compliquée à saisir. Fondamentalement, trois questions se posent à ce ni- 
veau : (1). En TQC, on exige que la dynamique de la théorie soit relativiste ; (2). Les va- 
riables de l'espace et du temps doivent être traitées aux mêmes pieds d'égalité, ce qui 
est précisément une autre exigence de la théorie de la relativité. En mécanique quan- 
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tique, le temps est juste un paramètre alors que la position est un opérateur (le symbole 
\t) est dénué de sens) ; (3). La mécanique quantique est avant tout une théorie à une 
particule, alors que dans la TQC plusieurs particules se mettent en jeu. L'intégrale de 
chemin répond admirablement à toutes ces exigences. 

Maintenant, voici le grand saut en un mot : la TQC remplace la coordonnée de po- 
sition x par un champ (p(x), où la quantité x est actuellement un raccourci pour les 
coordonnées x, y, z, t. Ce processus de réaffectation des coordonnées est connu sous le 
nom de seconde quantification. 

En TQC, nous supposons l'existence d'un champ quantique cp(x) qui pourra décrire 
des spécificités comme le spin des particules, moment angulaire, etc. Dans ce qu'on ap- 
pelle TQC canonique, le champ lui même est un opérateur (heureusement, les intégrales 
de chemin ignorent cette complication). A partir de ce que nous avons vu dans le cha- 
pitre précédent, il est tout à fait compréhensible que l'intégrale de chemin dans la TQC 
soit de la forme : 



où tous les coefficients libres sont actuellement "absorbés" dans la mesure *2)§ qui 
prend une forme comme d^(xj) à.tyix'i) • • • • Le fait d'appeler cette intégrale Z est une 
pure convention. D'une façon imagée, on pourrait voir cette intégrale comme l'ampli- 
tude de transition pour un champ qui se propage du "vide" à t = — oo au "vide" à 
t = oo. Un champ peut généralement décrire n'importe quelle dynamique, mais nous 
pouvons le voir dans notre situation comme une quantité décrivant une population 
de particules s'interagissant en chaque point de l'espace-temps. Il faut noter qu'appe- 
ler Z une intégrale de chemin n'est pas approprié puisqu'elle ne s'étend pas sur les 
chemins de l'espace-temps. On l'appelle l'intégrale fonctionnelle Z ou la fonctionnelle 
génératrice Z (generating functional). 

Parce que l'intégrale fonctionnelle avec termes d'interaction ne peut pas être évaluée 
directement, on a recours en général à une approche perturbative. 

La forme du Lagrangien pour un champ dépend de la nature des particules et des 
forces qui interviennent dans l'interaction. Par conséquent, il y a des Lagrangiens pour 
des particules scalaires (spin zéro) ou bosons, pour les spineurs (spin 1/2) ou fermions 




(4.1) 
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et pour des particules vectorielles (spin 1). Il y a même un pour les gravitons (spin 2). 
Le plus simple est le Lagrangien scalaire ou bosonique. Le présent chapitre s'illustrera 
sur une TQC scalaire relativiste ayant pour Lagrangien : 



L = -(du0 l <p-m 1 <p 2 )-V (4.2) 



- (dp*'' 

Le terme d'interaction le plus simple que l'on peut considérer est V ~ A^) 4 où A est 
une constante appelée la constante de couplage. C'est ce qu'on appelle une TQC à 
self-interaction (self-interacting field theory), autrement dit, le champ interagit avec 
lui-même et avec toute particule créée le long de l'interaction. Comme nous l'avons 
indiqué, l'intégrale Z n'a pas une forme explicite et de ce fait un traitement perturbatif 
s'impose. Ceci conduit à une interprétation fort intéressante du processus de création 
et de propagation des particules. 

Le problème est donc de résoudre l'intégrale 

Z(A) = /^expjj j - A(£ 4 d 4 xj (4.3) 

cette intégrale ne peut pas être attaquée dans cette forme. Le problème principal est que 
le nombre de variables d'intégration est infini. Nous aurons besoin de faire quelques 
changements avant de proposer une solution perturbative. 



4.2 Modification de la fonctionnelle génératrice Z 

Considérons la forme libre (A = 0) de Z avec une "source" J(x) : 




f tendra vers par la suite. Notons que nous nous plaçons dans un système d'unité où 
h = c = 1 et que dx remplace d 4 x pour alléger les écritures. Une intégration par partie 
du terme d^cpd^'cp donne 

J 1 (drfdt'cp - m 2 cp 2 ) dx = - J (^cpd 2 cp + \m 2 ^ 2 \ dx 
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Maintenant, supposons que le champ peut être mis sous la forme 

<p(x) = (p (x) + q>(x) 
où (po(x) est ce qu'on appelle la solution "classique" de l'équation hétérogène 

- (a 2 + m 2 )</>o(x) = / (4.5) 

c'est l'équation de Klein-Gordon avec source. La solution classique est unique (elle 
correspond au "chemin" classique). Elle peut être obtenue par la méthode usuelle des 
fonctions de Green : 

- (V + m 2 ) G(x-x') = ô i {x-x') (4.6) 

oùG(ï-x') est la fonction de Green (à 4 dimensions) associée à l'opérateur — (9 2 + m 2 ) . 
Donc 

(p (x) = j G(x-x')J(x')d A x' 

est la solution désirée. Pour résoudre (4.6) par rapport à G(x — x'), nous supposons 
que la fonction de Green peut être exprimée comme une transformée de Fourier d'une 
certaine fonction G(k) : 

G(x-x') = j -^G(k)e ik < x -^ (4.7) 

où k est le quadrivecteur énergie-impulsion ( E /c, ) . En faisant agir l'opérateur — (9 2 + 
sur l'expression de G(x — x'), et tenant compte du fait que la transformée de l'unité est 
tout simplement la distribution de Dirac : 

ôHx-x')= ! SLé k < x -*'\ 



(27T) 
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nous aboutissons à la solution 

G ( k ) = ÏJT^— 2' < 4 - 8 > 
k L — m L 

et par suite, 

G(x — x ) = («> 

(po(x) = J G(x-x')](x')A^x'. (4.10) 

Dans la TQC scalaire, la fonction de Green (4.9) est appelée le propagateur de Feyn- 
man et notée Ap(x — x'). Le propagateur G(k) exprimé dans l'espace des impulsions 
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nous sera utile plus tard quand nous définirons les règles dites de Feynman pour un 
processus de diffusion 1 . Pour récapituler, nous avons les formules : 



ou 



<Pq( x ) = j 


f A F ( 


x - x')J(x')d 4 x', 




Ap(x — x' 


> = J 


f d 4 Jc e ffc - (*-*') 


(2tt) 4 k 2 - m 2 ' 



Un calcul direct mais long permet d'écrire (4.4) sous la forme 

Z[J] = J@<p exp/^y d4j ~// K^OAf^-xOK^OdVdV 
Faisant sortir tout ce qui ne dépend pas de cp, on a 



Z[J] = Nexp 



J{x')A F {x" - x')J(x")d*x'd*x" } . 



ou 



N = J@cp exp \ l -J (d }l cpd }l cp - m 2 cp 2 ^ d 4 x 



peut être vue comme une constante de normalisation, mais nous allons poser N = 1 
parce que nous utiliserons des amplitudes normalisées plus tard. On a donc finalement 



Z[7] = ^p j -\ j/ }(x)A F (x - x')J(x')d 4 xd*x' 



(4.11) 



Nous avons donc transformé une intégrale à nombre infini de variables à une autre 
sur les quatre variables spatio-temporelles. Nous sommes maintenant en mesure de 
développer la fonctionnelle génératrice Z[J] pour un traitement perturbatif. 



4.3 Développement en série de puissances - les fonctions 
de Green 

La fonctionnelle Z peut se développer en série de Taylor comme suit (conformément à 
la formule (1.16)) 

Z t/] = E / - dx 1 dx 2 - ■ ■dx n J(x 1 )J(x 2 ) ■ ■ ■ J(x n )G(x lf . . . f x n ) (4.12) 

n=0 J U - 



1. Pour éviter la divergence de l'intégrale (4.9), le dénominateur de l'intégrand est généralement 
remplacé par k 2 — m 2 + ie. 
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ou 



G{%\, . . . , Xfi) — — 



Z[I] 



(4.13) 



_5]{x 1 ) SJ(x 2 ) à]{x n )\ " J /=o 
La fonction G est une sorte de fonction de Green. Dans la TQC elle est appelée fonction 
à n-points (n-point function). Nous allons voir que la fonction à n-points n'est non-nulle 
que pour des valeurs paires de n. Les physiciens ont compris que ces fonctions ne sont 
autres que les amplitudes de transition. Par exemple, dérivant (4.11), nous pouvons 
calculer la fonction à 2 points 



G(xi,x 2 ) 



Z-2 [ô}(x 1 )ë}(x 2 ] 

iA F (x 1 - x 2 ). 



z(7) 



/=o 



Schématiquement, Ap(x\ — x 2 ) peut se représenter par un graphe connexe à 2 points. 



^2 



Un calcul similaire donne pour la fonction à 4 points 

G(xi,X 2 ,X3,X4) = — Ap(xi — X 2 )Af(x3 — X4) — Af(xi — X3)Af {x 2 — X4) — Af(xi — X4)Af(x2 

ce qui peut se représenter par les trois graphes à 4 points 




Pour la fonction à 6 points, on trouve 15 termes correspondant à 15 graphes différents. 
En général, pour la fonction à 2n points, le nombre de graphes est le nombre de façons 
de grouper les 2n points en paires distinctes. Une analyse combinatoire inductive montre 
que ce nombre est égale à 



(2n - 1)!! = 1 x 3 x 5 • • • x (2n - 1) = ^ ^ = (4.14) 
K J v J 2 n - 1 (n-l)\ 2 n n\ 

Ce résultat peut se démontrer par récurrence. 



La fonctionnelle Z est un exponentiel. Le calcul peut donc être fait pour n'importe quel 
ordre. On verra qu'en général, le calcul est souvent laborieux. Heureusement, on verra 
aussi qu'il y a une formule simple pour ce genre de calcul perturbatif . 
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On interprète les fonctions G comme des quantités associées à des particules crées à 
un certain point de l'espace-temp X{, se propage le long d'une ligne (représentée par 
la quantité A]?(xf — xi)), et annihilée ensuite en un point final xt. Le propagateur de 
Feynman est le pivot des fonctions à n-points et de ce qu'on appelle les diagrammes 
de Feynman. 

Feynman a compris que les multiples differentiations (4.13) produit des termes que l'on 
peut décrire par des graphes simples. Chaque terme a un coefficient qui lui est associé 
relatif à la topologie du diagramme correspondant. L'ordre de multiplicité de chaque 
terme dépend du nombre de façons avec lesquelles on peut tracer le diagramme. Ces 
amplitudes et multiplicités peuvent être évalués par l'analyse combinatoire. 



4.4 L'intégrale fonctionnelle Z avec le terme d'interac- 
tion 



Nous allons maintenant introduire dans le Lagrangien le terme d'interaction A du po- 
tentiel (/> 4 . Par une convention qui sera claire par la suite, nous poserons V = A^ 4 /4!. 
La nouvelle fonctionnelle Z devient alors 



A 
4! 



d 4 x \ (4.15) 



Z(A) = [®tp exp jz j 1 (d }l (p{x)dV<p{x) - m 2 <p 2 (x 

= J@(p ex p{'/^ (d F (p(x)dt'(p(x) -m¥(x)j d 4 x-z^ J(p A (x)d 4 x 



Cette intégrale se présente vaguement sous la forme 



- uZ ~ bui du. 



Il est facile de montrer que l'on peut écrire cette intégrale sous la forme différentielle 
suivante : 



/ 



e~ u2 - bui du = exp 



Vc 



e~ u2+cu du 



c=0 



Nous rappelons que pour un opérateur O, exp(O) = Y^=q O n /nl. Nous allons user de 
cette approche pour introduire le terme exponentiel contenant A dans la fonctionnelle 
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Z[J]. Nous écrivons Z(A) comme 

Z(A) = exp 



.A y , 4 * 4 
*4! y X z^/(x)4 



Z[7] 



(4.16) 



Remplaçant Z[J] par son expression (4.11), nous obtenons 



Z(A) = exp 



.A r 4 ^ 4 
Z 4! y X ^/(x) 4 



exp l-^JJ J{x)A F (x - x')/(x')d 4 xd 4 x' 



(4.17) 



L'exponentiel de l'opérateur différentiel se développe comme suit 



exp 



.A r 4 Ô A 
" Z 4! ' d X ~ 



âj(x) 



1 - A (a* S * 
l — i—l a x 



4! J ô](xY 2!(4!) 2 



A2 V* si t* 5i 



E 

n=0 



-zA\" II- S A r Ô A 

— I dxi— — rj x / dx 2 —, — — x 



4! / n! 



SJ(x)*J ^/(y)4 
<5 4 



+ 



dx n 



(4.18) 



Étudions en détail les dérivées fonctionnelles de Z[J] par rapport à /. On a 

S 



ÔJ( Xl ) [ 
i 

~ ~2 



Donc 



l - j J(x')A F (x' - x")J(x")dx'dx" 

j ô(x' - xi)A F (x / - x")]{x")dx'dx" }{x')A F {x' - x")ô(x" - xi)dx'dx" 

= -i j A F (x 1 - x')J(x')dx' (4.19) 



ÔZ 



ÔJ( Xl ) 



■i J A f (xi - x')J(x')dx'^J Z. 



Nous avons de même 

Ô 2 Z 



-iA F (x 1 - x 2 )Z - (i J A F {xi - x')J(x')dx'J 



ÔZ 



SJ(x 2 ) 



SJ(x 1 )SJ(x 2 ) 

Il est commode d'adopter la notation concise suivante : 

Zx = SZ/SJfa), Z12 = ô 2 Z/ô]{x l )ô]{x 1 ),eïc A 12 = A F {x r - x 2 ) 
En éliminant l'intégrale, nous pouvons écrire 



Z12 = -z'A 12 Z + 2 Z i Z 2- 
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Seuls les termes qui sont proportionnels à Z survivent dans la limite / — >• 0. Il est aussi 
important de souligner que tous les termes avec nombres de différenciations impair 
s'annulent dans cette limite (comme Z12345). 

Pour évaluer la fonction à n points G(x\, . . . , x n ), nous devons évaluer des dérivations 
de Z(A) aux points X\, Xi,...,x n . Par exemple, la fonction à 4 points au premier ordre 
par rapport à A dans le développement (4.18) sera 



-xxxxl234 



/=« ÔJ(x)^J(x 1 )SJ(x 2 )SJ(x 3 )SJ(x 4 ) 



où x est le point d'interaction. La fonction à 4 points à l'ordre 2 est donnée par Z^^yy^^ | r 
Là, nous avons deux points d'interaction x et y, etc. 



4.5 Formules de diff érentiations et facteurs de symétrie 

Rappelons que la fonction à n points G{x\, . . . , x n ) se calcule par dérivation fonction- 
nelle de Z(A) n fois, évaluées pour / = 0. Pour la fonction à 2n points, nous avons déjà 
vu qu'elle se présente sous la forme générale 



- 123... 2n 



^ = {-i) n (A 12 A 34 • • • A 2n _i 2n + (Qi - 1) autres termes) (4.20) 



où le nombre C n est donné par (4.14). Comme nous l'avons expliqué, chaque terme 
peut être représenté par un graphe connexe. Ainsi, tout terme avec un nombre im- 
pair de dérivées fonctionnelles doit s'annuler vu que toute ligne dans le graphe doit 
joindre seulement deux points. La situation devient plus intéressante lorsque l'on doit 
considérer des propagateurs de la forme Ap(x — x) = Ap (0), où les points initial 
et final sont les mêmes. Par exemple, calculons Z xxxx correspondant à Z1234 lorsque 
x\ = xi = X3 = X4 = x. Dans ce cas, nous avons 



^ o = -3A f (0) 2 (4.21) 



Ryder [3] exprime graphiquement Ap(0) par un cercle Op our désigner qu'une parti- 
cule se crée au point x, se propage pendant un moment et ensuite se détruit au même 
point x, ce qui peut se présenter par une boucle ; l'analogie est très suggestive. Le 
terme (4.21) se schématise par -3QQ • Le facteur de pondération -3 est appelé le 
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facteur de symétrie pour Z xxxx . Sa valeur absolue représente le nombre de façons avec 
lesquelles un graphe peut être écrit. En général, quand le nombre de dérivations est 
assez grand, le calcul devient excessivement laborieux. Y-a-t-il un moyen d'écrire di- 
rectement toutes ces quantités ? La réponse est oui et c'est en plus relativement simple. 

De l'équation (4.20), nous voyons qu'en même point x, nous avons 

Z^_x = (-i) n C n A F (0) H . 

2n fois 

Le facteur de symétrie dans ce cas est donné par {—i) n C n . Mais qu'en-est-il lorsqu'on 
dérive en plus par rapport à des variables Xj qui ne sont pas nécessairement égales, par 
exemple 

Zxx\2 } Q = -A F (0)Ai2 - 2A xl A x2 
(où A x i = Ap(x — x\), etc). En général, dans la quantité 



Zxx . . . Xl2...2p 



2n fois 



(4.22) 



il y a, entre autres, des termes qui contiennent des facteurs Ap (0) k pour k < p ; les 
facteurs de symétrie correspondant à ces termes sont donnés par 

En fait, l'entier k est égale au nombre de boucles dans le graphe correspondant. Pour 
illustrer, un calcul de Z xxxx i234 1 / = o' ^ u ^' rappelons le, donne la fonction à 4 points 
évaluée à l'ordre 1 par rapport à A, donne 

^xxxxl23i 

= 3A f (0) 2 (A 12 A 34 + A 23 A 14 + A13A24) 

+ 12A P (0)(A rf A x2 A 3 4 + A x2 A x4 A 13 + A x3 A x4 A 12 

+ A x2 A x3 A u + A xl A xi A 23 + A xl A x3 A 2i ) + 24A. Tl A x2 A x3 A x4 . (4.24) 

Nous remarquons que les termes proportionnels à Ap (0) 2 , Ap (0) et Ap(0)° sont bien 
pondérés par (4.23) pour k = 2, 1 et respectivement, et de surcroît apparaissent avec 
des multiplicités respectivement égales à 3, 6 et 1. Ces multiplicités s'expliquent par 
le fait que les points x ; peuvent être permutés pour produire des graphes équivalents. 
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Les trois termes de la quantité (4.24) se schématisent par les diagrammes suivants 



X 2 Facteur de symétrie : 12 
x 4 Multiplicité : 6 






Facteur de symétrie : 3 
Multiplicité : 3 



Facteur de symétrie : 24 
Multiplicité : 1 



Une analyse combinatoire inductive montre que la multiplicité M pour tout terme 
dans (4.22) correspondant à un nombre de boucles k, s'écrit 

(2p)\ 



M 



in\(2p-2t)\ 



(4.25) 



ou 



t = k + p - n. (4.26) 
Il ne faut pas oublier que c'est une formule valable uniquement à l'ordre O(A). 

L'ordre n par rapport à A fait, en général, intervenir n dérivées quadruples par rap- 
port à X\, xj_,...x n . Ainsi, à l'ordre 2, on a à calculer Z xxxxlJl jy V i2...2p \ j =0 - Un exemple de 
graphes à 2 points, à l'ordre 2 est 
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La formule donnant le facteur de symétrie et la multiplicité pour un ordre n quel- 
conque (n vertexes internes) et 2p points externes, s'écrit 

n -i 

c = n 



W=l(i<7) ''■ 



(2p)\ 



2H\{2p-2t)\ 

n 

t = p-2n+ £ s i; -. (4.27) 
y=i(«'<;*) 

s,y désigne le nombre de lignes internes joignant les vertexes i et j. Par exemple, pour le 
diagramme précèdent, on a n = 2, p = 1, s X x — Syy = 0, s X y = 3, f = et M = 1. Donc 

c = I = I 

3! 6 



4.6 Règles de Feynman 

Pour calculer l'amplitude de transition pour un processus réel, Feynman a élaboré un 
ensemble de règles qui portent maintenant son nom. Rappelons (4.8), qui est le propa- 
gateur de Feynman dans l'espace des moments k. C'est effectivement beaucoup plus 
pratique que d'utiliser la forme spatiale A xy , et nous pouvons utiliser (4.8) pour asso- 
cier à chaque ligne dans un diagramme de Feynman avec un quadri-vecteur impulsion 
k. Les règles de Feynman pour la TQC <p^ sont 

1. Dessiner tous les diagrammes possibles correspondant au nombre désiré d'inter- 
actions (x, y, ...) et de points de l'espace-temps (x\, %i, ... ), utilisant le temps 
comme l'axe vertical et l'espace comme l'axe horizontal. Pour chaque graphique, 
étiqueter chaque ligne interne et externe avec une impulsion k et l'orienter par 
une flèche (le sens de la flèche peut être complètement arbitraire). Pour chaque 
ligne interne, écrire Y intégrale propagateur 

f dk i 
J {2nYk 1 -m 1 ' 

2. Pour chaque vertex d'interaction, écrire un facteur — ik. 
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3. Pour chaque vertex d'interaction, écrire une fonction delta de Dirac qui exprime 
la conservation d'impulsion sur ce sommet : la somme des 4-impulsions entrantes 
est égale à la somme des 4-impulsions sortantes : 



où k est positif si la ligne entre le vertex et négative si elle sort du sommet. 

4. Pour chaque graphique, il y aura une fonction delta résiduelle de la forme (2tc)^ô (k\ + 
kj_ + . . .) qui exprime la conservation globale de la 4-impulsion dans le diagramme. 



5. Dans l'évaluation des intégrales qui restent les fonctions delta introduites dans 
l'étape 3 simplifient énormément les calculs, et on n'a souvent même pas à les 
faire. 

6. Calculer le facteur de symétrie pour chaque graphique utilisant (4.27) et le mul- 
tiplier par le résultat obtenu à l'étape 5. 

Notons que pour tout processus physique, il peut y avoir un très grand nombre de 
schémas possibles en fonction du nombre d'interaction envisagés. Quand l'ordre de 
l'interaction augmente, le nombre et la complexité des schémas possibles augmente 
rapidement. Toutefois, la constante d'interaction est généralement une quantité très 
faible (dans l'électrodynamique quantique, il est numériquement égale à environ 1 / 137), 
la petitesse de ce terme pour n grand réduit efficacement la probabilité qu'un processus 
compliqué se produise effectivement. Cela justifie l'utilisation de l'approche perturba- 



Exemples. Pour voir comment ces règles fonctionnent, nous allons voir deux exemples 
Pour le diagramme 




Éliminer ce terme. 



tif. 




x 



4 



X 
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Ici n = 1 et, par conséquent, C = 1 et M = 1, l'amplitude est donc simplement —z'A. 
Maintenant, calculons l'amplitude pour le diagramme à l'ordre 2 suivant 




Nous avons deux points d'interaction qui contribuent avec le facteur (— z'A) . Nous 
avons deux lignes internes, donc nous avons 

z = HA)2 / \ ^^k^^h^^ 5 ^ + k z- k -ri w 4 ^ k + « - *3 - m 

simplifiant, on trouve 

Z = A 2 / 1 7 Tl r^-Tyy .(27T) 4 ^ 4 (fc 1 + k 2 - k 3 - fc 4 ) 

i (2zr) 4 k 2 - m 2 (ki + k 2 - k) 2 - m 2V ' v y 

Il y une fonction de Dirac résiduelle qui traduit la conservation de la 4-impulsion total. 
On peut s'en débarrasser. Pour ce diagramme, s xx = s xx = 0, s xy = 2. Utilisant (4.27), 
on a 




L'amplitude est donc 

1 2 f d 4 fc 1 1 

2 7 (2tt) 4 (fc 2 - m 2 ) (fci + fc 2 - fc)2 - m 2 

Mais, est-ce que cette intégrale converge ? NON ! ! !, car pour k — > oo, l'intégrand se 
comporte comme logfc — >• oo. C'est le célèbre problème de la divergence ultraviolette. 
La solution repose sur une acrobatie qui est vivement critiqué pour son manque de "ri- 
gueur mathématique" selon les mots de Paul Dirac. Cette technique marche pourtant 
bien et donne des résultats satisfaisant ; elle est connue sous le nom de renormalisation. 
Malheureusement le temps réservé à ce cours ne permet pas de s'étaler sur ce point. 
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